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Abstract. The anomalously large rates of some hadronic transitions from quarkonium are studied using QCD multipole 
expansion (QCDME) in the framework of a constituent quark model which has been successful in describing hadronic 
phenomenology. The hybrid intermediate states needed in the QCDME method are calculated in a natural extension of our 
constituent quark model based on the Quark Confining String (QCS) scheme. Some of the anomalies are explained due to the 
presence of an hybrid state with a mass near the mass of the decaying resonance whereas other are justified by the presence 
of molecular components in the wave function. Some unexpected results are pointed out. 
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INTRODUCTION 

The recent experimental data of the hadronic transi¬ 
tions of heavy quarkonia, such as y/(n5) or Y{nS) to 
lower states with emission of two pions, show a puz¬ 
zling behavior. The experimental data show that the 
A (4260), A(4360) and A(4660) states apparently de¬ 
cay only through particular channels. The A(4260) res¬ 
onance has been seen only in J/y/n^n^, whereas the 
A(4360) and A(4660) appear only in the \j/{2S)7t^7r^ 
channel. Furthermore, these two resonances have an 
anomalously large width [1, 2, 3]. Moreover, such decays 
have not been observed for the ij/(4415). In the bottom 
sector, compared to the ordinary T {nS) —> Y (mS) {m <n) 
transitions, the partial width of the T( 10860) is out of 
line by two orders of magnitude [4]. 

The typical momentum involved in the transition is too 
low to use perturbative QCD and therefore such meth¬ 
ods can not be applied. Then, non perturbative methods, 
like the QCD multipole expansion approach (QCDME), 
should be used to describe at least the transition of 
the lower lying states. In the single channel picture of 
QCDME the heavy quarkonium system serves as a com¬ 
pact color source which emits two soft gluons which 
hadronize into two pions. After the emission of the first 
gluon and before the emission of the second one there 
exists an intermediate state where the QQ pair together 
with the gluon forms an hybrid state. The width of the 
transition critically depends on the particular spectrum 
of the hybrid states, therefore it is important to describe 
the QQ states an the hybrid consistently using as few pa¬ 
rameters as possible. 

Apart from lattice calculations [8, 9], hybrid me¬ 
son properties have been calculated in different models: 


the flux-tube model [10, 11], constituent gluons [12], 
Coulomb gauge QCD [13] and quark confining string 
model (QCS) [14, 15, 16] or QCD string model [17]. 
Among them we adopt the QCS model since it was used 
in the early works of QCDME and it incorporates finite 
quark mass corrections. 

Above the open flavor threshold, the QCDME scenario 
may change due to the possible contribution of molecular 
components [5], which can modified the hadronic decay 
width. 

In this work we will address the description of the new 
data of the hadronic transitions in heavy quarkonium us¬ 
ing the QCDME in the framework of a constituent quark 
model (see references [6] and [7] for reviews) which has 
been successful in describing the hadron phenomenol¬ 
ogy and the hadronic reactions. Hybrid states are con¬ 
sistently generated in the original quark model using the 
QCS scheme. In this way we minimize the number of 
free parameters. Above the open flavor threshold molec¬ 
ular meson-meson channels coupled to the cc states are 
included in the calculation 


THEORETICAL FRAMEWORK 

The constituent quark model we use is based on the as¬ 
sumption that the spontaneous chiral symmetry break¬ 
ing generates the constituent quark mass. To compensate 
this mass term in the Hamiltonian the Lagrangian must 
include Goldstone boson fields which mediates the inter¬ 
action between quarks. The minimal realization of this 
mechanism include a pseudoscalar boson and an scalar 
one. This fact does not affect the heavy quark sector but 
is of paramount importance in the molecular picture be- 


cause the only remaining interaction between the two 
molecular component, due to its color singlet nature, is 
the one driven by the Goldstone boson exchanges be¬ 
tween the light quarks. Below the chiral symmetry scale 
quarks still interact through one gluon exchange and con¬ 
finement potential. Explicit expression for this interac¬ 
tions are given in Ref. [19]. 

In this model mesons are described as clusters of 
one pair of quark and antiquark. To found the quark- 
antiquark bound states we solve the Schrodinger equa¬ 
tion using the Gaussian Expansion Method [20]. In the 
case of the molecular structures the four body problem 
is also solved using the gaussian expansion and the two 
body wave functions obtained from the solution of the 
Schrodinger equation. To derive the meson-meson inter¬ 
action from the qq interaction we use the Resonating 
Group Method (RGM) [21]. The coupling between the 
qq and the four quarks configuration are performed us¬ 
ing the model [22]. A more detailed discussion of the 
model and its application to the states mentioned above 
can be found in Ref. [5]. 

To calculate the hadronic transition we use the 
QCDME method. The multipole expansion has been 
widely used for studying radiation processes in which 
the electromagnetic field is radiated from local sources. 
In our case the energy difference between the initial and 
the final states is usually small and therefore the gluon 
wavelength is large compared to the typical size of the 
QQ states and can be treated in a multipole expansion. 
The gauge-invariant formulation of multipole expansion 
within QCD was given by T.-M. Yan in Ref. [23]. 

The two pion hadronic transitions are dominated by 
double electric-dipole transitions (ElEl). The transition 
amplitude split into two factors. The first one concerns 
to the wave functions and energies of the initial and fi¬ 
nal quarkonium state as well as those of the intermediate 
states. All these quantities can be calculated using suit¬ 
able quark models. The second one describes the con¬ 
version of the emitted gluons into light hadrons. Its scale 
is the light hadron mass, which is very low and there¬ 
fore cannot be calculated in a perturbative way. Usu¬ 
ally one uses a phenomenological approach based on 
PCAC [24] which involves two parameters (Ci and C 2 ). 
The Cl term is isotropic (S-wave) while the C 2 term is 
angular dependent (D-wave). These two parameters are 
fitted to the well established ^f(2S) J and 
V/(3770) J transitions. 

As explained above, these transitions involve the emis¬ 
sion of two gluons that hadronized into two pions. The 
intermediate states after the emission of the first gluon 
and before the emission of the second gluon are states 
with a gluon and a QQ pair which are the so-called hy¬ 
brid states. We describe these states in the QCS model. 

The QCS model is defined by a relativistic-, gauge- 
and reparametrization-invariant action describing quarks 


interacting with color SU (3) gauge fields in a two dimen¬ 
sional world sheet. The model has no gluonic degrees of 
freedom, but has instead string degrees of freedom. 

The string can carry energy-momentum only in the 
region between the quark and the antiquark, thus the 
quarks appear to be at the ends of the string. The equa¬ 
tion that describes the dynamics of the quark-antiquark 
pair linked by the string is the usual Schrodinger equa¬ 
tion with a confinement potential. Gluon excitation ef¬ 
fects are described by the vibration of the string. These 
vibrational modes provide new states beyond the naive 
meson picture. A complete description of the QCS model 
can be found in Ref. [15]. Its generalization to our con¬ 
stituent quark model is described in Ref. [25]. 

RESULTS 

All the parameters of the quark model are taken from 
Ref. [26]. The quantum numbers of hybrid states which 
participate in the two pion transition are = 1 . 

The mass of the two low lying states in the ccg sector 
are 4.35 and 4.64GeV. The first one agrees with the 
results of the flux-tube model [11] (4.1 —4.2), Coulomb 
gauge QCD [13] (4.47), QCD string model [17] (4.40), 
potential models [18] (4.23) and lattice calculation [9] 
(4.40). In the bottom sector the ground state mass is 
10.785 GeV near to the mass of the potential model 
(10.79) [18] 

Table 1 shows the calculated ^jt+n-y/{nS) ^ ^e+e- 
for the = 1 charmonium states. As the decays 
V/(25) J/xj/n+n^ and ^r{2>110) J have 
been used to fit the Ci and C 2 parameters, they are not 
included in the table. One can see that in the case of the 
decay channel xi/{2S)n^n^ the only significant values 
correspond to the decays of the X (4360) and X (4660) 
which are also in agreement with the recent experimen¬ 
tal data. 

In the decay channel Jjxj/n^n^ a high value of the 
^jc+n-y/{nS) ^ U+f- obtained for the Y(4360) reso¬ 
nance. This result apparently contradicts the experimen¬ 
tal data because this decay has not be reported in the reac¬ 
tion > y/ xj/n^n^ [1]. The cross section of this re¬ 

action shows a resonance in the 4.2 — 4.4 GeV. energy re¬ 
gion which has been attributed to the X (4260). This res¬ 
onance does not appear in our calculation as a cc meson 
and its nature is still under discussion. An interference 
between the X (4260) and X (4360) resonances would be 
possible and should be explored. 

If one looks to the values of the cross section at peak 
(Table 2) one sees that in the decay channel J/xj/n^n^ 
the values of the rest of peak cross sections, besides the 
resonance at 4.36 GeV, are of the order of the experi¬ 
mental background and it is difficult to decide whether a 



TABLE 1. Ryi^ns) = ^n+n-\ir[nS) X '^e+e- for the = 1 charmo- 
nium states. Experimental data are from Ref. [3]. 


Initial Meson 

R‘h^(eV) R5)^(eV) 

Ki2S) ^Vi2S) 

Vr(4040) 

1.20 

0.11 

Vr(4160) 

0.40 

6 X 10^2 

A (4360) 

52.5 

5.05 7.4 ±0.9 

A(4415) 

3 X 10^® 

0,27 

A (4660) 

0.58 

1.08 1.04 ±0.5 


TABLE 2. The cross section at peak for the = 1 S-wave char- 
monium states. Experimental data are from Ref. [3]. 


Initial Meson 


^y7v'(pb) 

<(25) (Pb) 

<(25) (Pb) 

Vr(4040) 

13.46 

- 

1.25 

- 

Vr(4160) 

3.32 

- 

0.50 

- 

A(4360) 

329.7 

- 

31.69 

52±2 

A(4415) 

4 X 10-5 

- 

3.38 

- 

A(4660) 

10.97 

- 

20.29 

28±2 


TABLE 3. R-Y(nS) — ^K+n-i:{nS) X ^e+e- for the = 1 bottomonium states. Experimental data are from Ref. [27]. 


Initial Meson 

<(!,) (eV) 


<(2S) 

RTiis) 

<(3S) 

<(35) (eV) 

T(25) 

r(35') 

T(45) 

T( 10860) 

98.34 

23.94 

6 X 10-2 

4.1 X 10-2 

105.4±4.3 
18.5±9.8 
(2.3 ±0.9) X 10-2 
1.64 ±0.40 

5.58 

2.5 X 10-2 
5.8 X 10-2 

(2.3 ±0.4) X 10-2 
2.42 ±0.64 

1.8 X 10-2 

1.49 ±0.65 


resonance is present or not. These results justify why no 
signal of the v/(4040), V/(4160) and v/(4415) has been 
reported. 

The results for the bottomonium sector are shown in 
Table 3. One can see that the theoretical values agree 
reasonably well with the experimental ones except in the 
case of the T( 10860). In this case we do not find any 
hybrid state near its energy and the mechanism which 
explains the large width in the charm sector cannot be 
applied to this case. 

The T(10860) is above the open-bottom threshold and 
then contributions of molecular components of the wave 
function like BB, BB* or B*B* may appear. The impor¬ 
tance of these components lies in the fact that the de¬ 
cays of the molecular structures through the Ynn chan¬ 
nel is OZI allowed whereas the decay of a bh component 
is OZI forbidden. Therefore the presence of molecular 
components in the wave function enhances the transition 
probability. 

Following Ref. [5], we have performed a coupled 
channel calculation including the bb 5^Si and A^D\ to¬ 
gether with the BB, BB* and B*B* molecules. The re¬ 
sult of this calculation is shown in Table 4. We obtain 


two resonances in the 10860 MeV mass region. The first 
one is too narrow and the second is wide. We iden¬ 
tify the wide resonance with the T( 10860). The process 
T(10860) —>■ Y{nS)nn involves the creation of a qq pair 
and a quark rearrangement process. This can result into a 
nonresonant final state T {nS)nn, or an intermediate state 
in which the two pion forms a resonance (like /o(980)) 
following by a subsequent decay into the final channel. 
In the calculation we include the /o(500) and the /o(980) 
resonances together with the nonresonant contribution. 
Results are shown in Tables 5 and 6. The contribution of 
the molecular components is able to explain the branch¬ 
ing ratio of the T(10860) —> T(15);r;r but fails in the 
Y(2S)nn and Y{3S)nn cases. In all cases the main con¬ 
tribution corresponds to the nonresonant process but this 
contribution decreases as the phase space decreases. The 
contribution of more resonances could help to improve 
the results. 

SUMMARY 

Using the constituent quark model quoted above we cal¬ 
culate the J/\j/nn and v/(25);r;r decays ofthe/^*'= 1 



TABLE 4. Molecular components of the 1 bottomonium states above the 
BB threshold.. 


Mass (MeV) 

Width (MeV). 

Prob (5^5i) 

Prob (4^Di) 

Prob(Mol) 

10829.10 

75.31 

0.72 

11.27 

88.01 

10849.73 

15.51 

0.24 

2.38 

97.38 

10849.93 

89.43 

3.48 

11.99 

84.53 


TABLE 5. Widths and Branching Ratios of the Y{nS) 
hadronic transitions through molecular components. Ex¬ 
perimental data are from Ref. [27]. 


Channel Width (MeV). B.R.[%] B.R. exp [%] 

T(15);r;r 0.330 0.362 0.53 ±0.06 

T(25);r;r 0.023 0.026 0.78 ±0.13 

T(35);r;r 0.004 0.004 0.48 ±0.19 


TABLE 6. Resonant and Nonresonant contribution to the 
T(nS) hadronic transition . 


Channel /o(500) /o(960) Nonresonant 

T(15);r;r 5.27 x lO^^* 1.14x10^3 0.330 

T(25);r;r 9.22 x lO^^* 1.30 x lO^^ po3 x lO^^ 

T(35);r;r 5.97 x 10^"^ 2.86x10^^ 1.98 xlO^"* 


charmonium and bottomonium states. Hybrid states are 
consistently generated in the original quark model using 
the quark confining string (QCS) scheme. Above thresh¬ 
old we incorporate the effects of molecular components. 
We are able to explain the anomalously large decay width 
of the A(4360) and A(4660) due to the presence of hy¬ 
brid states located near the masses of both resonances. 
However this mechanism does not work in the case of the 
T( 10860). The contribution of molecular components 
improves the description of the T(10860) -5- Y{lS)nn 
branching ratio but fails to explain the branching ratio in 
the Y{2S)nn and Y{3S)nn cases. Our model also pre¬ 
dicts a large value for the decay A (4360) J/xj/nn 
which has been not seen in the experiments. Therefore 
the puzzling situation of the hadronic decays still de¬ 
serves more experimental and theoretical efforts. 
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